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We apply optimized perturbation theory to the quark-meson model at finite temperature T and 
quark chemical potential /i. The effective potential is calculated to one loop both in the chiral limit 
and at the physical point and used to study the chiral dynamics of two-fiavor QCD. The critical 
temperature and the order of the phase transition depends heavily on whether or not one includes 
the bosonic and fermionic vacuum fluctuations in the effective potential. A full one-loop calculation 
in the chiral limit predicts a first-order transition for all values of yU. At the physical point, one finds 
a crossover in the whole /i — T plane. 
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I. INTRODUCTION 



■ If we consider QCD in the chiral limit, the QCD Lagrangian has an extended symmetry. For Nf quark flavors, the 
II ' Lagrangian is invariant under SU (A^/)l x SU {Nf)ii transformations. In the vacuum, this symmetry is spontaneously 
I- broken down to SU{Nf)v by the formation of a chiral condensate. According to Goldstone's theorem, this gives 
^ ' rise to Nj — 1 massless modes, the Nanibu-Goldstone bosons. 

In the real world, quarks are not massless. If we add a mass term for the quarks to the QCD Lagrangian, 
chiral symmetry is explicitly broken, which means that the Goldstone bosons become massive as well. Nevertheless, 
the u and d quarks are light enough that we have an approximate chiral symmetry in Nature, and for simplicity, 
K>" ' we therefore often consider them to be massless. QCD with two flavors of massless quarks is a commonly used 
approximation, the validity of which seems to be supported by the observed lightness of the pions. However, there 
are more than two quark flavors in the real world. In particular, the strange quark mass is on the order of 100 McV, 
and so we expect it to have an impact on physics at this energy scale. 
• , If one heats strongly interacting matter, it is expected that there will be a transition to a phase where chiral 
' symmetry is restored and the quarks within the hadrons are liberated. The location and nature of the transition 
^— ^ ' from the phase with broken chiral symmetry to the chirally symmetric phase is a topic of active research. The order 
I of the phase transition depends on the number of quark flavors and their masses. For example, for two quark flavors 
• • ] and vanishing baryon number chemical potential fi — 0, we expect, a second-order phase transition in the chiral limit 
' if the U{1)a is explicitly broken, otherwise it is driven flrst order. In the case of flnite quark masses it is a crossover. 
For three massless quarks, Nf — 3, the transition is alwaw first order. These results follow from universality-class 
5^ , arguments and the (non) existence of stable fixed points 1, 2]. At nonzero baryon chemical potential and in particular 
at T = 0, the order of the chiral transition is not obvious from universality arguments [3j. However, most model 
calculations predict a first-order transition at T = 0, 3 ■ If this is the case, there is a line of first order chiral 
transitions in the /i-T plane that ends in a second-order transition. 

Currently, the only way to obtain quantitative information about the chiral phase transition from first principles 
is to use lattice Monte Carlo simulations. Results from lattice QCD indicate that the transition takes place at 
Tc ^ 150 MeV depending on the number and masses of the quarks 0, 0. 

Unfortunately, lattice QCD is severely limited by the sign problem. When the baryon number chemical potential 
is nonzero, the fermion determinant, which is used as the statistical weight in the Monte Carlo importance sampling, 
becomes complex. There are several methods to sidestep the sign problem for small fi. Examples include reweighting 
techniques, Taylor expansions in fi around /z = 0, and the use of imaginary chemical potentials. For a review, see 
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Ref. [7}. However, these methods only work for fJ. ^T, and a general solution to the sign problem for all ^ has yet 
to be found. 

A complementary approach to the investigation of the chiral phase transition, which works even at nonzero baryon 
number chemical potential, is offered by effective models — simplified models that share QCD's chiral symm etry 
breaking pattern. Examples include linear sigma models [Ml, NJL-tiT)e models [iMl- Polyakov-NJL models^ 
the quark-meson (QM) effective model ^ [la. [l9l-[28j . (see also Ref. Q), and the Polyakov-QM model [Ullc" 
In this paper we will employ the QM model. 

By definition, when using effective models we neglect certain aspects of the full theory. In xPT and sigma models, 
for instance, we only have mesonic degrees of freedom. For this reason, there is no way to include a nonzero baryon 
number chemical potential in these models. 

In NJL-type models we have quark degrees of freedom, however, the color symmetry is global, which means that 
there are no gluon degrees of freedom. One-gluon exchange has been replaced by local four-fermion interactions 
between the quarks. The NJL model is not confining and the quarks are not asymptotically free 4]. Unlike the 
scalar models above, it allows the inclusion of separate chemical potentials for different quark fiavors (and colors), 
and through the use of NJL models a rich structure of color superconducting phases sd, fi ^ and low temperature 
has emerged (ssj . 

While the effective models themselves are already simplifications of a more complicated theory, one often makes 
additional simplifications within the models. For qualitative considerations, for instance, one may decide to omit 
contributions that do not seem immediately important for the problem at hand. As an example, in the QM model, 
chiral symmetry breaking takes place in the meson sector, and as a consequence, the vacuum contribution to the 
free energy from the fermions is sometimes omitted. In the NJL model, on the other hand, this term is responsible 
for the chiral symmetry breaking, and so it cannot be neglected. 

When making such simplifications, however, it is important to know exactly what one is discarding, in order 
not to "throw the baby out with the bath water". To wit, for the example just given, it was recently shown that 



neglecting the fermion vacuum contribution to the QM free energy changes the order of the phase transition [26 1. 
In the present paper we will study the role of the vacuum fiuctuations in detail and generalize to /i 7^ 0. 

The quark-meson effective model contains, as the name implies, both quark and meson degrees of freedom. The 
mesons are described as in the linear sigma model, by an 0(7V)-symmetric scalar field theory with a quartic self- 
interaction. The meson sector couples to the quark sector through Yukawa-type couplings. 

We shall investigate the thermodynamics of the QM model using optimized perturbation theory (OPT), which 
was introduced in 1998 by Chiku and Hatsuda @, In this framework, a local quadratic term proportional 

to an arbitrary parameter x is added to and subtracted from the Lagrangian d, [13, l40l - |48| In the perturbative 
expansion, the added mass term is included in the free part of the Lagrangian and the subtracted term is treated 
as an interaction term. This constitutes a reorganisation of the perturbative series in which higher-order terms of 
nai've perturbation theory are resummed to all orders. 

We have already mentioned that the pions are the Nambu-Goldstone bosons associated with the spontaneous 
breaking of chiral symmetry. As such, the Nambu-Goldstone theorem plays a significant role in the study of the 
chiral phase transition. Importantly, OPT ensures that this theorem is satisified order-by-order of the expansion 

0- . . 

The paper is organized as follows: In Sec.|TTl we present the quark-meson effective model and the OPT framework; 
in Sec. IIII[ we calculate the free energy to one loop; in Sec. IIV[ we present our results and in Sec. El we make a few 
concluding remarks. Finally, expressions for various sum-integrals are provided in an appendix. 



II. QUARK-MESON MODEL 

We consider the quark-meson effective model with N real scalar fields ((/)i, (j)2, ■ ■ ■ , i^n) and Nf massless quark 
flavors {ip) with Nc colors. The bare Euclidean Lagrangian for the model can be written in the form 

Cb^C + AC, (1) 



^ The QM model is also known as the linear sigma model coupled to quarks (LSMq). 
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where £ is the renormalized Lagrangian and A£ contains the counterterms. For clarity, we further write £ as 



£ — -ClSM + £q + /^Yukawa- (2) 

The first term in this expression, £lsm, is the Lagrangian of the hnear sigma modeh 

ChSM = ^{dti4>i)idf,(j)i) + ^ml(l)i(j)i + - h(j)i . (3) 

In the case of vanishing /i, this Lagrangian has an 0{N) symmetry which, assuming TOq < 0, is spontaneously 
broken down to 0{N — 1). When h ^ the symmetry is explicitly broken, mimicking the explicit breaking of chiral 
symmetry in the QCD Lagrangian by nonzero quark masses. 
Next, we have £q, describing the massless quarks: 

£q = i/'(7m^m - M74)V' ■ (4) 

Here, fi is the quark number chemical potential, which we assume to be equal for all quark flavors, and which is 
related to the baryon number chemical potential hy fi — ij.b/5. Cq has a U{1)b symmetry related to the conservation 
of baryon number, as well as an SU{Nf)L x SU{Nf)i^ chiral symmetry. 

In Eq. (ID) and throughout the remainder of the paper, we take 7^ to denote Euclidean gamma matrices. These 
are related to the familiar gamma matrices in Minkowski space by 7,; = i"f\,j and 74 = j^.j. They have the following 
properties: {7,1,71/} = 26^,^1 and tr[7^7y] = AS^^. 

^Yukawa dcscribcs the Yukawa couplings between the scalars and the quarks, and we write it compactly in the 
form 

^Yukawa = gi^Ti^cjii. (5) 

Here, we have defined 

r=(l,-*75r), (6) 

where r = (ti, T2, T3) are the Pauli matrices. This only makes sense when = 4, which is the value we will use in the 
end. The fifth Euclidean gamma matrix is defined as 75 = 71727374 = ~7m- Finally, the counterterm Lagrangian 
A£ is written as 

1 1 CA 

A£ = -A{^M^^.'I^^) + l^Bmlcj),cj,, + -^{M^f - Dm^ + ^^7^^^^ + G#r,V'</>, , (7) 

where A, B, C, D, F, and G are the renormalisation constants for the scalar field, the scalar mass, the four-scalar 
coupling, the vacuum energy, the fermion field, and the Yukawa coupling, respectively. In the present calculation, 
we need to know B, C, and D to leading order in the couplings. 



A. Broken symmetry 

The group 0(4) is homomorphic to SU{2) x SU{2). The consequence of this is that if we take = 4 and Nf = 2, 
an 0(4) rotation among the scalar fields can always be countered by an SU{2) x SU{2) transformation on the 
fermions, keeping the entire Lagrangian invariant. 

Furthermore, 0(3), the symmetry group of >Clsm in the broken phase, is homomorphic to SU{2), which is the 
symmetry group of Cq after chiral symmetry has been broken by, say, the addition of a nonzero mass term. Indeed, 
we shall see that spontaneous symmetry breaking in the mesonic sector in fact induces a dynamic mass term for the 
fermions. 

In the phase with broken chiral symmetry, the scalar field (j) takes on a nonzero vacuum expectation value, and 
provided h > 0, the vacuum state will always point in the (f)i direction. When h = 0, 0{N) symmetry allows us to 
choose any direction for the vacuum state, and for simplicity we choose the same direction in this case. We therefore 
write 

(j) = {a + v,n), (8) 
where a and tt are fluctuations representing the sigma meson and the pions, while v is the vacuum expectation value 

of 
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B. Optimized perturbation theory 



In optimized perturbation theory (OPT), we add and subtract quadratic terms in the Lagrangian: 

n 1 2 1 1 2 1 

L + -XCr + -XTTiTTi - -X(T - -XTTiTTi . (9) 

The first two terms are treated as part of the free Lagrangian and make this expUcit by absorbing them into a new 
mass parameter m, defined by 

= jtIq + X ■ (10) 

The subtracted terms are treated as two-particle interactions. Optimized perturbation theory is then defined by the 
counting rule that an insertion of x counts as a factor of A or in loop diagrams. 
Note that we must take care to adjust the counterterms as well. For instance, 

A£o = ^B{m^ - xV + - D{m^ - xf ■ (11) 

The optimization parameter x is in principle completely arbitrary, and if wc could carry out calciilations to all 
orders, the results would be independent of x- In order to complete a calculation, we therefore need a prescription for 
X- There are infinitely many possible choices, but we limit ourselves to mentioning the following two prescriptions: 

• The principle of minimal sensitivity (PMS) says that one should try to minimize the dependence of the result 
on the parameter x- Thus, a PMS condition on some observable calculated to i loops, Og,, would then be 
expressed as 

• A fastest apparent convergence (FAC) condition means that the perturbative corrections in 0(, should be as 
small as possible. Sometimes it is even possible to enforce the strict condition that 

Oi - Ol-n = (13) 

for some n satisfying 1 <n < I. 



III. FREE ENERGY 



We now calculate the free energy to one loop. In the following, we employ the notation 




(15) 

for the dimensionally regularized sum-integral, where d = 3 — 2e is the dimension of space and A is an arbitrary 
momentum scale. The factor e^ /A-k is chosen so that, after subtraction of the poles in e, A coincides with the MS 
renormalization scale. 

To make the various contributions to the prefactors of the diagrams as explicit as possible, we take the number 
of scalar fields to be N (i.e. (TV — 1) pion fields) and the number of quarks to be Nf. 
At tree-level, the free energy is given by 

^0 = ^^m'v^ + j^y^ - hv . (16) 
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To one-loop order, we obtain 



Tl = J'la + -Fib + -Flc - ^X^^ + A J"l 



where 



and AJ^i is the one-loop counterterm: 
The masses are 



i^log(p2-fm2), 



2 I 2 \ 



-NcNfJl tr log(«7pP^ + mq - ^74) 



(17) 

(18) 
(19) 

(20) 
(21) 

(22) 
(23) 

- 5-«- ■ (24) 

The relevant sum-integrals are given in Appendix [X] Expanding the temperature-independent terms to order e*^, 
we find 



2,^2 
m + —V 



2,^2 
m H V 

6 



,2„,2 



^Ib 



1 



A 



/ A 

(47r)2 \mq 



1 3 

1 3 

7+2 



+2Jo(/?m,)r4 
-I- 2Mfim^)T^ 



1 3 

7+2 



where Jo{(3m) and Kq {l3mq,(3fj,) are defined in Appendix A. The divergent terms in J^ia-ic are thus 



1 



la— Ic 



(25) 
(26) 
(27) 

(28) 



4(47r)2e 

In the minimal subtraction scheme the counterterm AJ^i cancels exactly the divergences in the free energy (|28 
which yields the one-loop counterterms 

iV + 2 A 



6(47r)2 e 
1 

(47r)2e 
N 1 



i(iV + 8)A-24AreA^/^ 



(29) 
(30) 
(31) 



4(47r)2 e' 

The final result for the renormalized one-loop free energy is obtained by adding Eqs. (^5]) - (|27p and Eq. (I?!]) using 
Eq. dMll-dSI]) 



(iv-i),4- 



47VeiV/m^ 



+2 [Jo(/3m,) + {N~ l)Jo(/3m^) + 2A^,A^/i^o+(/3m„ /S^i) -f 2N,NfK^ {iSnig, T' 



(32) 
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where we have defined 

La = log a = a,TT,q. (33) 

ml 

The full one- loop free energy J^o+i is then given by the sum of Eqs. p6|) and (p2|) . Using the renormalization group 
equation for the running coupling A, one can show that J^o+i is renormalization group invariant. 



IV. RESULTS AND DISCUSSION 



In this section we discuss the procedure for performing numerical calculations using OPT, and we present our 
results. 

At one-loop order, the PMS condition on the free energy, 



, (34) 



dx 

simply has no solution. 

Instead, we consider the gap equation for v, obtained by differentiating the free energy: 



dv 

Here, Hi is the one-loop pion self-energy at zero external momentum 



v{mi + ni)-h = 0. (35) 



nr = i^. (36) 

V dv 



A FAG condition on the self-energy is then Hi = 0, or 



2 



where 



Sim, T, ^i) = [L„, + 1] + -^[K+{/3m, ^3^i) + K^i(3m, /3^i)] 



(38) 



This condition on x has a straightforward physical interpretation — it means that Goldstone's theorem is fulfilled 
even at tree level: 

ml = -. (39) 

V 

This also means that we can write = h/v + and thus all ^-dependence is removed from the r.h.s. of Eq. 

1571 in the broken phase. 

In order to perform numerical calculations, we need to find numerical values for the model parameters mg. A, g, 
and A in terms of the vacuum {T = fi = 0) values of m„, m^r, mq, and v. Matching at tree level yields 

^l + X = -\{ml-?.ml) , (40) 

<) , (41) 

(42) 
(43) 



A = 


3. 2 
V 




iriq 


9 = 






V 


h = 


m^v. 



^ In some cases 111 = has no solution. We then minimized |ni| instead in accordance with the FAC condition. 
^ At tree level, it is the combination m? = 4- X that is determined, cf. Eqs. II22I I- II23II . 
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At tree level, A is undetermined. We choose the renormalization scale to be 

n2 



to: 



Finally, we set iV = 4, Nt = 2, and = 3. 



A. Contributions to J^o+i 

At the one-loop level we distinguish between the following contributions to the free energy: 
1. The tree-level contribution from the bosons, 



x-trcc -77 



-m?v^ + — Tj"* — hv. 
2 4! 



2. The vacuum contribution from the bosons, 

1 



" ~ 4(47r)2 
3. The thermal contribution from the bosons, 

■J h 



{N - l)mi 



rp4 



2(47r)2 

4. The vacuum contribution from the fermions, 



-77vac _ ^ J ' 

f - (47r)2 



5. The thermal contribution from the fermions. 



J" 



(47r)2 



[K+{^im,,(i^J) + K^{(3mq,(3^l)\ 



(44) 



(45) 



(46) 



(47) 



(48) 



(There is also the OPT term, — i which we omit for the purpose of the following discussion. It is taken into 
account in all numerical calculations.) We will now investigate the effect of these various contributions on the 
structure of the phase diagram. 



Chiral limit 



We start by considering the phase diagram in the chiral limit, h — where Lagrangian has an exact chiral 
symmetry and the pions are massless Nambu-Goldstone bosons. We then fix our model parameters using the 
vacuum values 



and using Eqs. (gOll-dlll), we find that 



= 600 MeV , 

TO^ = , 

WLq = 300 MeV , 
V = = 93 MeV , 

+ X = -(424.26 MeV)2 

A = 124.87, 

g = 3.2258, 

/i = , 

A = 363.92 MeV . 



(49) 
(50) 
(51) 
(52) 

(53) 
(54) 
(55) 
(56) 
(57) 
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Starting with the simplest case, we may say that we are only interested in thermal effects, and thus ignore the 
vacuum contributions J^™'^ and J'J^'^- Furthermore, it is common to use a mean- field approximation to the meson 
sector, in which we ignore the thermal fluctuations, J"^, as well. Hence, all that remains of the free energy is the 
tree-level boson contribution and the thermal fluctuations of the quarks: 



/ 



(58) 



The phase diagram one obtains using this approximation is shown in Fig. [T] [left panel). We see that the phase 
transition is first order for all values of /i, a fact which is explicitly demonstrated for = in the right panel. The 
critical temperature is Tc — 140 MeV, which is reasonably close to the critical temperature of QCD. The critical 
chemical potential for T — is /ic — 296 MeV. 




100 150 200 250 300 
Quark chemical potential [MeV] 




120 130 140 

Temperature [MeV] 



FIG. 1: At one loop, including only the thermal contribution of the fermions and the tree-level contribution of the bosons: 
(left panel) the phase diagram in the — T plane, where the solid line represents a first-order phase transition, and (right 
panel) the solutions of the gap equation plotted as a function of T a,t fi — 0. 



Recently, Skokov et al [26| showed analytically for /i = that if one includes the fermion vacuum term, i.e. 

V , (59) 



the order of the phase transition changes — it becomes a second-order phase transition. For T = and finite /i, this 
can be shown analytically as well. This is displayed in the right panel of Fig. [2l Furthermore, as the left panel 
shows, the phase transition is of second order for all values of ^. The critical temperature at = is Tc = 190 MeV, 
a significant increase as compared to the case where the fermion vacuum fluctuations were ignored. Similarly, there 
is a large increase in the critical chemical potential for T = Q to jic = 345 MeV. 




100 150 200 250 300 
Quark chemical potential [MeV] 




1 70 1 80 1 90 
Temperature [MeV] 



FIG. 2: At one loop, including both the vacuum contribution and the thermal contribution of the fermions, but only the 
tree-level contribution of the bosons: (left panel) the phase diagram in the fJ^ — T plane, where the dashed line represents a 
second-order phase transition, and (right panel) the solutions of the gap equation plotted as a function of T at /i = 0. 
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Next, let us omit the fermionic vacuum fluctuations again, and instead study the effects of including the thermal 
fluctuations of the mesons: 



Trtrcc 
^ h 



(60) 



This situation was investigated by Bowman and Kapusta in Ref. j25| . In the chiral limit, they found that there is a 
first-order phase transition extending all the way from the fi axis to the temperature axis, with a critical temperature 
of Tc ~ 140 MeV at = 0. We show this in Fig. [31 where we also see that we get a lower critical temperature of 
about 110 MeV. The critical chemical potential is for T = is fic = 357 MeV. 




50 100 150 200 250 300 350 
Quark chemical potential [MeV] 




1 00 110 
Temperature [MeV] 



FIG. 3: At one loop, including the thermal contributions from both fermions and bosons, but not vacuum contributions: 
(left panel) the phase diagram in the fi — T plane, where the solid line represents a first-order phase transition, and (right 
panel) the solutions of the gap equation plotted as a function of T at = 0. 

We have now seen that by including J^™'^ the first-order phase transition turns into the second-order transition 
one would expect at high temperatures, and that by taking J^^ into account we pushed the critical temperature 
down to the vicinity of T^'~^^. It is therefore tempting to guess that by including all contributions. 



(61) 



we will get a second-order phase transition near T = 150 MeV. Unfortunately, this turns out not to be the case. As 
Fig. U shows, the critical temperature is 148 MeV, but the phase transition is still first order — in fact it is even more 
strongly so. Part of the reason for this is the fact that the vacuum energy of the mesons has opposite sign from the 
vacuum energy of the quarks and so the contributions cancel to some extent .The critical chemical for potential for 
T = is 357 MeV. 




100 150 200 250 300 
Quark chemical potential [MeV] 




135 140 145 
Temperature [MeV] 



FIG. 4: At one loop, including all contributions from both fermions and bosons: (le]t panel) the phase diagram in the \i. — T 
plane, where the solid line represents a first-order phase transition, and (right panel) the solutions of the gap equation plotted 
as a function of T at = 0. 
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C. Physical point 



In the chiral limit, v is exactly zero in the chirally restored phase. If we let h > 0, chiral symmetry is explicitly 
broken in the Lagrangian. Then, the chiral condensate never vanishes completely, but approaches zero asymptotically 
at high temperature and/or large baryon chemical potential. Furthermore, the pions acquire mass and become pseudo 
Nambu-Goldstone bosons. 

At the physical point, we use the vacuum values 



nicr 

m„ 



and using Eqs. pOjl - P^ we find the parameters 

ml + x 
A 

9 
h 
A 



600 MeV , 
139 MeV , 
300 MeV , 
U = 93 MeV 



-(388.611 MeV)2 
118.17 , 
3.2258 , 

(121.573 MeYf , 
363.92 MeV . 



(62) 
(63) 
(64) 
(65) 



(66) 
(67) 
(68) 
(69) 
(70) 



Let us look at what happens to the phase transition in the case when only thermal contrubtions to the free energy 
are taken into account. 




50 100 150 200 250 300 350 
Quark chemical potential [MeV] 



FIG. 5: At one loop, including all contributions from both fermions and bosons (upper curve) and thermal contributions 
only (lower curve). 



This is shown in Fig [S] (lower curve). The line of first-order transition that starts at fic — 287 MeV and T ~ ends 
at a critical point at ^ — 220 MeV and T — 90 MeV. In comparison, Scavenius 1131 et al. obtain /j, = 207 MeV and 
T — 99 MeV for the position of the critical endpoint, while Bowman and Kapusta f25| obtain the values /i — 283 MeV 
and T = 75 MeV. 

If we include the bosonic as well as the fermionic vacuum contributions, i.e. we include all contributions in the 
model, the critical endpoint disappears. The transition is a crossover in the whole /i — T plane. This is shown in 
Fig [S] as well for comparison (upper curve). Only for lower values of m^r will the critical endpoint appear. 



V. SUMMARY AND OUTLOOK 



In the present paper, we have applied optimized perturbation theory to the quark-meson model to study the chiral 
transition in the fi — T plane. The properties of the phase transition - the order and the critical temperature - depend 
crucially on the the approximations made. The one-loop calculation in the chiral limit yields either a first-order or a 
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second-order transition in the whole ^ — T plane depending on whether one includes the fermion vacuum fluctuation 
term or not. This is in contrast to most NJL-type model calculations Q and functional renormalization g roup 
calculations that predict a curve of first-order transitions starting at T = ending at at critical point [2l|, [l^, l36l |. 
Taking a model seriously means including the effects of all its degrees of freedom. It therefore questionable to throw 
out some terms unless one can show that they are not important. For example, the effective potential is receiving 
contributions from vacuum fluctuations from all energies scales up to the ultraviolet cutoff A of the effective theory. 
This is the region of validity of the low-energy effective theory and it is therefore not obvious that they can be 
omitted. In view of this, it is perhaps somewhat surprising that the critical endpoint disappears all together at the 
physical point and we are left with a crossover in the whole ^ — T plane. 

We close with a few remarks concerning future work. Optimized perturbation theory was used at the two-loop 
level in (/)'*-theory to study the phase transition and it was shown that it correctly reproduces the second-order nature 
at zero baryon chemical potential [10|]. It would therefore be of interest to extend the two- loop calculations to the 
quark- meson model at finite /i along the lines of Ref. [l^. Secondly, one could apply OPT to the case Nf — 3, where 
we know the transition in the chiral limit is first order. Dialing the strange quark mass from = cx) to ttt-s = 
interpolate between the chiral limit for Nf — 2 and Nf — 3. Somewhere along the way, the transition then changes 
from second order to first order and it would be interesting to find the critical value for the strange quark mass, m* 
where this happens. 



Appendix A: One-loop sum-integrals 

We need the following bosonic sum-integrals 

^log {P'^+m') 



1 /A 



(47r)2 \m 
1 /A 



2e 



2e 



-.(l-.)(2-e)"^ -Jo(/3™)T 



The temperature-dependent term is expressed as integrals over the Bose-Einstein distribution function: 



Jnix) 



4eTTfi) 



oo i4-2ri-2£ 



r(l-n-e) JO 
These integrals satisfy the recursion relation 

dx X 
We need the following fermionic sum-integrals: 



dt 



1 



2xJn+iix) 



(Al) 
(A2) 

(A3) 

(A4) 



V tr log(i7^P^ + m- ii'ji) 
T{P} 



1 /A 



(47r)2 \m 



pTerCl _|_ f] 



Here, the temperature-dependent term is expressed as integrals over the Fermi-Dirac distribution function: 



oo j.4-2n-2e 



dt 



1 



(A5) 



(A6) 



The functions K^{x,y) satisfy the same recursion relation as Jn{x), 
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